Thermal transport in one-dimensional spin heterostructures 
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We study heat transport in a one- dimensional inhomogeneous quantum spin 1/2 system. It 
consists of a finite-size XX spin chain coupled at its ends to semi-infinite XX and XY chains at 
different temperatures, which play the role of heat and spin reservoirs. After using the Jordan- 
Wigner transformation we map the original spin Hamiltonian into a fermionic Hamiltonian, which 
contains normal and pairing terms. We find the expressions for the heat currents and solve the 
problem with a non-equilibrium Green's function formalism. We analyze the behavior of the heat 
currents as functions of the model parameters. When finite magnetic fields are applied at the two 
reservoirs, the system exhibits rectifying effects in the heat flow. 
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I. INTRODUCTION 

In the last decade there has been a renewed inter- 
est related to the research of thermal transport in one- 
dimensional magnetic systems^ 2 - . Most of these studies 
have been motivated by unusual high values of thermal 
conductance of some materials, as for example reported 
in Ref. H. From the theoretical side, several works cal- 
culating and discussing thermal conductivity in different 
model Hamiltonians have also appeared™. On the other 
hand, spin Hamiltonians provide the natural scenario to 
implement quantum computation. This motivated in- 
teresting proposals of a variety of physical systems, like 
arrays of quantum dots, optical lattices and nuclear mag- 
netic resonace experiments which are architectured to ef- 
fectively behave like one-dimensional spin systems. 6 

Most of the theoretical studies on heat transport in 
spin systems are performed within the linear response 
regime, assuming a very small temperature gradient X/T, 
and using a Kubo formula. Altough this formula is widely 
used in calculating thermal transport properties it has 
several conceptual difficulties, particularly if \jT is not 
small-i The thermal averaging assumes a constant tem- 
perature, and (as in all Kubo formulas), the time depen- 
dence is governed by the total Hamiltonian H = Hq + H', 
including the perturbation H' ~ i while the ther- 
modynamic averaging is done with Hq, assuming a fast 
evolution. For finite the separation of H into an 
unperturbed part Hg and a perturbation H' is in princi- 
ple ambiguous. On the other hand, real experiments are 
performed by coupling the system under study to macro- 
scopic systems with well defined temperatures that act 
as reservoirs. Therefore, it seems desirable to develop 
alternative approaches designed to treat systems out of 
equilibrium. Important progress has been achieved by 
studying the properties of non-equilibrium steady states 
of XX and XY chains within an algebraic setting which 
allows one to obtain explicit analytical expressions for dif- 
ferent quantities such as entropy production^. Recently, 
the usefulness of Kubo formula for the investigation of 
heat transport in quantum systems has been discussed^ 



and the investigation of energy transport through spin 
systems beyond Kubo formula has been addressed on 
the basis of master equations and quantum Monte Carlo 
methods^ Nevertheless, it is not clear how to extend 
these ideas to more general situations. 

In this work we address the problem of thermal trans- 
port in one dimensional magnetic systems from a differ- 
ent perspective. We study a problem that to the best of 
our knowledge has not been considered yet. We study the 
heat flow through a spin chain heterostructure, which we 
generically define as a set of finite or semi-infinite spin 
chains attached by their ends. Each piece of the het- 
erostructure can be in principle described by a different 
Hamiltonian. For dcfiniteness, we shall consider in this 
work a finite central system connected to two (left and 
right) semi-infinite chains, with a finite temperature dif- 
ference between the m 10 ! 11 . This problem is the thermal 
analog of electronic transport through mesoscopic struc- 
tures, nanodevices or molecules connected to conducting 
leads with a finite applied bias voltage, a subject of in- 
tense research in recent years. In fact, this type of set-up 
is the common situation found in the study of charge 
transport in electronic systems, where a central system 
is connected to charge reservoirs that also act as thermal 
baths. This is the basis of the "Landauer" approach, 
which is one of the most common frameworks to study 
transport properties of meso- and nano- devices in the 
last years ji 2 - 

In addition, as it is well known, one-dimensional spin 
1/2 systems can be mapped, via the Jordan- Wigner 
transformation to fermionic systems. Thus, the model 
under investigation is equivalent to an electronic het- 
erostructure where very well established techniques, as 
the Schwinger-Keldysh non equilibrium Green's func- 
tions method can be applied^. 

We focus on a simple device where the central and right 
parts of the system can be described by XX spin 1/2 
chains, and the left part corresponds to an anisotropic 
XY chain. The Jordan- Wigner transformation maps 
these models into bilinear fermionic systems, rendering 
the theoretical study simpler. We show that this simple 
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FIG. 1: (Color online) Schematic digram showing the set 
up of the system. The right and central chains correspond to 
XX models while the left chain is described by an XY model. 
Left and right chains are taken semi-infinite at temperatures 
Tl and Tr respectively. 



device presents an interesting physical effect: due to a 
mechanism reminiscent of Andreev reflections in super- 
conductors, this device could act as thermal diode. This 
rectifying effect might be useful for applications. We ar- 
gue that this is a generic feature, which remains valid for 
anisotropic XYZ models. 

The paper is organized as follows: in Section 11 we 
present the model, and the non-equilibrium formalism 
based on Green's functions. Section 111 contains the nu- 
merical results. Section IV is a summary and discussion. 



II. THE FORMALISM 
A. Model 

We consider a system of three one-dimensional spin 
1/2 chains coupled through their ends as shown in Fig. 
[TJ The system is described by the following Hamiltonian 
H = Hl + He + Hn 



H C oup, where 
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where in the first line, the index a = L,C, R (left, center, 
right) labels the different chains. Each chain is character- 
ized by its nearest- neighbor exchange J a , with anisotropy 
A a for the ratio of the interaction along the y direction 
with respect to that along x, and the magnetic field B a 
applied along the z direction. The isotropic case for a 
given chain corresponds to A Q = 1. The Hamiltonian 
H C oup describes the coupling between the central chain, 



and the left and right ones. In terms of the representa- 
tion for the spin operators: S ± — S x ± iS v , the above 
Hamiltonian reads 
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For the isotropic case, X a = 1 only the flip-flop terms 
with products of one raising and and one lowering spin 
operators survive. 

We now introduce the Jordan- Wigner transformation 
to map the spin 1/2 Hamiltonian into a fermionic Hamil- 
tonian through: S+ t = ftjexp(inJ2'j fl,jf<*,j), where 
denotes a sum over all the positions located at the 
left of the position a, I. Similarly, the other spin oper- 
ators transform as = exp(-wrX^- fl,jfa,j)fa,l and 
S a,i = fljfa.l - V 2 : where the operators //,// obey 
fermionic commutation rules: {fa,i,fa>i'} = ^l,l'^a,a', 
and {fl v f a ,y} = {f a:h f a ', i'} = 0. Substituting in the 
Hamiltonian |2]), we get 

N a -i 

H a = ^ i w <x[f a ,lfcx,l+l+fi > l +1 f<x,l\ 

1=1 

+^aifa,lfa,l+l + fa,l+lfa,l]} 



A*a E] ^cn , I fa , I ) 
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[fllfc.l + fcjL.l] + 

w '[fc,N fR,i +/b,i/c,/v c ], 



(3) 



where 



= Ja(l + A a )/4, A Q = J a (l - Aa)/4, 
li a = B^, and w' = J'/2. Therefore, in the language 
of fermionic operators, the Hamiltonian contains "nor- 
mal" terms, with a creation and a destruction operator, as 
well as "anomalous" terms, with two creation or two de- 
struction operators. The normal ones are a hopping term 
between nearest neighbors (w a ), which is originated in 
the flip-flop spin terms and a chemical potential (/i a ) cou- 
pled to the fermionic density, which comes from the mag- 
netic field pointing along the z-direction. The anomalous 
terms (A a ) are similar to those of a one-dimensional 
Hamiltonian with a gap function with p-wave symme- 
try, decoupled in the Bardeen-Cooper-Schrieffer (BCS) 
approximation and are originated by the anisotropy be- 
tween the X and Y exchange interaction. 
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Inspired in this analogy, we focus our study on a junc- 
tion between a chain with isotropic interactions (XX 
spin chain) and an anisotropic one (XY spin chain), 
which in the fcrmionic language is similar to a normal- 
superconductor junction. Such a situation is realized in 
a configuration with \r = Ac = and Al ^ 0. We 
also assume that the left and right chains are at tem- 
peratures Tl and Tr respectively and they are both of 
infinite length( Nl — > oo and Nr — > oo) 

B. Energy balance 

A consistent procedure to define an expression for the 
heat current from first principles, is to analyze the evo- 
lution of the energy stored in a small volume of the sys- 
tem and derive the corresponding equation for the con- 
servation of the energyJ^ For the present Hamiltonian 
we choose an elementary volume containing two nearest- 
neighbor positions of the chain. We place the volume 
enclosing the sites 1,1 + 1 within the central (XX) chain, 
which in the fcrmionic language contains only normal 
terms. We work in units where ft = 1. The equation for 
the conservation of the energy enclosed by this volume is 

dEi : i +1 _ Jc d + + d(Sf) 

= w c-^(fl,ifc,i+i + fc,i+ifc,i) - 

= -iw c {[H, fljc.i+i + fc,i + Jc,i\) 

+i»c{[Hjy Cl i}) = J? +u+2 - J ; Q u,(4) 

where </^ +1 is the heat current flowing from I to I + 1, 
which in the present setup coincides with the energy cur- 
rent. Its explicit expression is obtained from the evalua- 
tion of the above commutator, which gives 

• J — J i,i+i — J i-iu 

= i£ u+i( £ ?-i,i(fcj-ifc\i+i - fc,i+ifc,i-i) 
+£?+i,i+i(fc,i+2fc.i+i - fc,i+ifc,i+2)), (5) 

where ef v denotes the matrix element Z, /' of the Hamil- 
tonian He- I n order to evaluate the above current it is 
convenient to introduce the lesser Green's functions 

G^AtJ)=i(fl v (t')Mt)), (6) 

thus 

J Q = 2Re{sf l+1 Gc l+lcl _ 1 (t,t)efL ll 

+ £ i+l,i+lCci + l,Ci+2(^' t) £ ?+2d+l\- CO 

The lesser Green's functions are one of the basic ele- 
ments within Keldysh non-equilibrium Green's function 



formalism-^. They are evaluated by solving the equations 
of motion (Dyson's equations), which for our model can 
be written as follows 

E G ^c fe M[4,^- £ ^,] = o (8) 

k 

for coordinates I, I' lying within the central chain. We 
have used the stationary property of the system, as 
a consequence of which the Green's functions depend 
on the difference t — t', which allows us to transform: 
G tA f - = I^dco/(2n)e-^ t ">G< 3 ,(cj). Thus, us- 
ing the above equation in Eq. ([7]) the energy current can 
be also expressed in the following way 

/ + OO 7 
— ^ +1 G< +liCi M 
-oo z7r 

+ £ ?+id+i G ci+i,ci+i( UJ )]}- ( 9 ) 

However Re{G^ +lj+1 (o;)} = 0. Thus, the heat current 
reduces to 

/-(- OO 7 
^ ^ £ £ +1 G< +1 ^M}. (10) 

Similarly, if we evaluate the heat current through the 
contacts L — C and C — R, we find 

/+oo J 
-oo 27r 

/+OC J 
-oo 27T 

C. Solving Dyson's equations 

In order to evaluate G < and the heat current we follow 
a treatment close to that presented in Ref. ITU. We define 
the retarded "normal" and "Gorkov" Green's functions 

F%At,f) = -m-mft^flvV)})- ( 12 ) 

Before writing down the Dyson's equations satisfied 
the these "full" Green's function let us define the following 
"free" particle g°(w) and hole 3 a (w) Green's functions 

CM]i=^("+il) + ^. (13) 

with 7/ = + . 

From now on we will work in Fourier space and we will 
not write explicitly the u> dependence of Green's func- 
tions unless necessary. For the left chain, we also define 

the functions G° L , G L containing the paring term contri- 
bution, through the relations 

[Gir 1 = [gir 1 - A L f L a l , 

P L ]- l = \f L ]- l -A L f L k L , (14) 
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Fl = G A L g u L , 
f L = G°A L f L . 



(15) 



Here ef^ and A a ; al , are matrices defined on the coordi- 
nates of the chain a = L, C, or R, containing respectively, 
the normal and anomalous elements of the Hamiltonian. 
In the case we are studying only is non- vanishing. 

To obtain these Green's functions, we write the fol- 
lowing Dyson's equation which relate them with Green's 

functions of the "disconnected" chains <?° and ~g a (see 
below) and the matrix elements of the contacts 



Mr 1 + i^r 1 + [<&] 

[W}}G R - A L F R = 1, 



{[fir 1 H 

[W]}F R 



fa}' 1 
-A L G { 



[9r 

= 6. 



(16) 
(17) 



The matrix W — Wl + Wr contains the matrix elements 
of H cont describing the connections between the central 
and left parts and between the central and right parts. 

The above equations can be rewritten in a more conve- 
nient form by recourse to the following procedure. From 
Eq. (HZD 



F R =f[A L G I 



WF ± 



(18) 



G R = g° + g°[A L F R + WG R ], (19) 
Substituting JTHJ in HH) and (fT9|) in (fTTf one obtains 



G R = G°(l + WG H ) + F L WF l 



F R = Fl(l + WG") -GWF 



R.\ 7=7 



?R 



(20) 
(21) 

Let us now consider Eq. (|20|) for the following partic- 
ular coordinates 

^CltCV — 9c,l,l' ~r 9c,l,l w u Ll,Ct' 

+9h,* w ' G RL,ci" ( 22 ) 

r R f^O nnT R 

U L1,CI' — Lz L,l,l u; U C1,CI' 



-Fl^w'F^ 



CVi 



U R1,CI' — 9r,1,1 w ^CN.CV 



(23) 
(24) 



Substituting Eqs. (USD and (J21J) in Eq. $22$ it is easy to 
see that the Dyson's equation for the two indices corre- 
sponding to coordinates of C can be written as follows 

{[gO]- 1 - £ R ^}G R + t R ^F R = 1, (25) 

{[f"]" 1 - £ R >ff}F R + t R >^G R = 6, (26) 



where the matrices of the above equations have sizes 
Nq x Nc and elements corresponding to the coordinates 
of the central chain. The "self-energy" matrices are 

Eg/' = S u ,\w'nS L {G° Lxl + S liNc g% hl 

Eg," = & iI VI 2 Pi,iG£ 1 i,i + w&,i j i] 

Eg/* = StMwf^F" (27) 



The explicit expressions for these functions imply the 
evaluation of all the functions appearing in the right hand 
sides of (|27| . Notice that these functions have been de- 
fined from manipulations of the Dyson's equations cor- 
responding to Hl or H r isolated from the central chain. 
We indicate a procedure for the calculation of these func- 
tions in appendix A. Note also that since the right and 
left parts of the system are held at two different but 
constant temperatures, these Green's functions can be 
calculated at equilibrium. The advantage of the above 
representation becomes clear by writing (f26|) as 



F R =g- c £ R ^G R , 

Iter 1 = [fcT 1 



(28) 



and substituting it in (|25|) . The result leads to the solu- 
tion of the retarded normal Green's function within C 



g r c = mr'-t 



J eff 



: R,9f± c £R,f a 



(29) 



The results obtained so far correspond to the retarded 
Green's functions and self energies. The lesser Green's 
function with coordinates within C can be easily ob- 
tained by recourse to Langreth rule o 16 i 17 . In particular 
one obtains^ 



c- 



(30) 



where the advanced Green's function is obtained from the 
retarded one, by means of the relation G c {uj) = [G^(u)Y 
and the lesser component of the "self-energy" is 



£<=£<'" +A R £ <j9 
eff 

±<^A A + A R £<<"A A , 



(31) 



with A r {uj) = £^Ht c (w) and A a (uj) = [A^(w)]t. 
The self-energies have components 

E^f'" (w) = i5i,i'[5i,\T v £ v (uj)f L (uj) +8i, n ^ V r (w)/r(w)] 

being T v £ u '(u) = -2Jm\Z R ' vv ' and T R y {oj) = 
-2Im[S R ' vu ' {lu)n,n} with v,v' = g,f. The Fermi func- 
tions f a (uj), with a = L, R depend on the tempera- 
tures Tjj and Tr of the left and right chains respectively: 
f a (w) = 1/(1 + e w / T °), in units where fcs = 1. 
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Finally, the lesser counterparts of Eqs. (|23f and (|24[). 
which correspond to Green's functions with one of the 
coordinates in the central (C) chain and the other one 
in the left (L) or right (R) chain, can be calculated by 
recourse again to Langreth rule s 16 ! 17 

Gli,ci' = ^l%x w ' ^ci,ci' + Gl^,i w 'Gci,ci< 

+* f> £,i«/Fci,cv + ^l^'Fcx,cv > (32) 

G< 0,< I f^iA , 0,-R / s~i< fQQ\ 

ri,cv ~ 9R t i t i w u cn,ci' ~t~ 9r,i,i w u cn,ci'- v-J-J; 

D. Heat currents and transmission functions 

We focus on the expression for the heat current eval- 
uated in the contact between the central chain C and R 
given in Eq. ([Til . Using Eq. ([55]) one obtains 

/+OO 7 

+Gg W)CW (w) fl S< 1 (w)]. (34) 
Using (j3TJ)) and after some algebra (see Ref. [HI), it is 



found 






JQ = 


/•+OO J 

/ — - / A (w)][T"(a;) 


- T a (u)}, 


T n (w) = 


rfM|Gg^>)| 2 r^u,), 




r» = 


rffMlA^H^r^H, 


(35) 


where 








rf + 2Re[rl^ A ,] + |A^ 1 | 2 r{^ 


" \ a i,n\ 2t r> 


A N,N — 




(36) 



The difference of Fermi functions in the expression of 
J® , reflects the fact that the existence of a non- vanishing 
heat current through the central system depends on the 
existence of a difference of temperatures between the left 
and right chains. The details of the model are enclosed 
in the behavior of the "normal" and "anomalous" trans- 
mission functions T n {ui) and T a (u), which are analogous 
to those defined in Ref. [TH in the context of particle 
transport in a setup with normal and superconducting 
wires. The first function has, in fact, the structure of a 
transmission. Notice that it depends on the densities of 
states of the right and left chains through the functions 
Tr and T e /f , and one the particle propagator between 
the first and last points of the central chain. Instead, 
T a (cu) actually has the structure of a reflection process. 
Notice that it depends on the density of states for par- 
ticles and holes of the right reservoir and on a multi- 
particle propagator N at the last point of the central 
chain. Typical plots for these functions are shown in Fig. 
[2J These functions do not depend on the temperatures 
Tl and Tr and are non-vanishing only within a finite 
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FIG. 2: (Color online) "normal" (left) and "anomalous" 
(right) transmission functions. The latter is the counterpart 
to Andreev transmission function in electronic systems. Pa- 
rameters are Al = 0.2, Jl = Jr = Jc = J' = 1 and all 
chemical potentials set to zero. 



range of energies of a width that is set by the largest 
exchange parameter between the left, right and central 
chains. These functions are symmetric with respect to 
u> = for Bl = Br = (see Fig. [2]). This symmetry is 
broken for finite B a , since the effect of a finite magnetic 
field in one of the side chains is to shift the corresponding 
function as T™' (w) -► T™' (u - B a ). 

In the language of fermionic systems, two different 
kinds of processes take place in a normal-superconductor 
junction. For energies higher than the gap, the trans- 
port is due to the tunneling of normal single particle 
high-energy excitations. This mechanism contributes to 
the electronic transmission function T n (u>). Instead, for 
low energies, below the gap, the transport is due to 
the mechanism known as "Andreev reflection" , which 
implies the combination of two fermions of the nor- 
mal side into a Cooper pair within the superconducting 
one, leaving a hole that is reflected back from the junc- 
tion into the normal side. Because of this mechanism, 
T n {ui) ~ T a {uj) ~ 1 for energies within the supercon- 
ducting gap, i.e. \u)\ < Al. The effective conversion 
of electrons into Cooper pairs taking place in the mech- 
anism of Andreev reflection helps to partice transport. 
Mathematically, this is reflected by the fact that the to- 
tal particle transmission function is T n (ui) + T a (ti>)j 15 i 18 
Instead, in the case of heat transport, T a (uu) and T n (uj) 
contribute with opposite sign, as explicitly shown in Eq. 
(f3"5)l . i.e. the mechanism of Andreev reflection, plays a 
negative role regarding the heat transport. The conse- 
quence is a vanishing heat transport due to excitations 
within the energy window defined by the superconduct- 
ing gap. 

In the original language of interacting spins, the above 



picture translates as follows. Low energy spin excitations 
traveling from the isotropic chain via flip-flop processes 
in the z direction meet an energy gap at the other side of 
junction due to the anisotropic interaction which tends 
to favor flip-flop processes in a different direction. This 
favors the simultaneous raising or lowering of two spins 
at two neighboring positions of the chain and causes mul- 
tiscattering processes in which a portion of the incident 
spin wave packet manages to twist and propagate into 
the other side, at the same time that a portion becomes 
reflected and propagates back. 

We can describe the behavior of J® for low T and small 
temperature gradients ST as follows. Writing Tr = T 
and Tl = Tr + ST we can approximate the difference of 
Fermi functions in Eq. (|35| as 



(37) 



/lM - f R {oj) gf~ 

On the other hand, from Fig. [21 we can write, 



r( w ) -T a H 

T n {uj) -T a {uj) 



leading to, 



ST 

7T 



uj < 2A L 

1 2A L < uo < 2 J 



2Al 



dT 



(38) 
(39) 



(40) 



For low enough temperatures, T -C 2A^ <C 2 J, this ex- 
pression can be further approximated as, 



k oT J 2Al 



-ST 

7T T 



Aj _2A £ 



(41) 
(42) 



Therefore, for T < A/,, the heat current is exponen- 
tially small. 

On the other hand, for A L = 0, the behavior of the 
is fully due to normal tunneling. For low T we can 
perform a Sommerfeld expansion on the Fermi function 
to get 



39 - 

~ ^TST, T«J Q , A i= 0. 



(43) 



III. RESULTS 



In this section we discuss the behavior of the heat cur- 
rent as a function of the different ingredients of the spin 
system. For simplicity, we consider identical exchange 
parameters along the left, central and right chains: 3 L = 
Jc = Jr = J' = J ■ Without loss of generality we set 
J = 1. Thus, we focus on a spin heterostructure with a 




FIG. 3: (Color online) Top left: Conductance, defined as the 
ratio of heat current to temperature gradient between right 
and left chains. Top right: specific heat of the central system. 
Bottom: ratio of the above magnitudes. Parameters are Tl = 
T + 8t,Tr = T where St = 0.005, all chemical potentials set 
to zero, and (from top to bottom) Al ~ 0,0.25, 0.5. 



single junction between a semi-infinite XX and a semi- 
infinite XY chain, which in the fermionic language trans- 
lates to a single S-N junction. For this particular config- 
uration, our results do not depend on the length of the 
central chain. 

As discussed in the previous subsection, the structure 
of the expression ([33)) for the heat current clearly reflects 
the fact that for small temperature differences, we obtain 
a behavior of the form 



J Q = -K5T, 



(44) 



where the coefficient K can be interpreted as a thermal 
conductance. It is tempting to relate this coefficient with 
the conductivity k evaluated in several works on the basis 
of linear response theory. If we assume that the relation 
between the two coefficients is similar to the one between 
electrical conductance and electrical conductivity, K and 
k should differ just by a geometrical factor. However, to 
the best of our knowledge, a rigorous relation between 
these two coefficients has not been presented in the liter- 
ature. Nevertheless, the behavior of if as a function of 
T shown in the left upper panel of Fig. [3J for the case of 
two connected XX chains (see the plot corresponding to 
Al = 0) is similar to the one reported in the literature for 
homogeneous and isotropic chains^^. In this case, the 
anomalous component is zero and K increases linearly in 
T for low temperature [see Eq. (|43p], as discussed at the 
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FIG. 4: (Color online) Heat current as a function of the 
temperature of the left-chain (right chain at fixed Tr = 0) for 
all chemical potentials set to zero, all hopping parameters set 
to 1 and different values of the paring parameter At, from 
top to bottom At = 0.1, 0.25, 0.5, 0.75. 

end of the previous section. The conductance reaches a 
maximum at T ~ J. and decreases at higher temper- 
atures, as a consequence of the finite bandwith (energy 
window) for the spin excitations amenable to cross the 
central chain transporting energy from one side to the 
other one. As expected, for a fixed temperature K de- 
creases for increasing values of At- In agreement with 
the behavior discussed in the previous section, K is ex- 
ponentially small for T < Al [see Eq. (l40l) ]. For higher 
temperatures, the high energy excitations are allowed to 
perform tunneling above the energy gap, with the con- 
comitant increase of K. As in the case with At = 0, the 
maximum is achieved at T ~ J. 

We also evaluate the specific heat for the equilibrium 
central system in contact to the side chains at the same 
temperature T as follows: 

i=i • / -°° 

In a normal metallic system as described by Drudc 
model, this quantity is related to the thermal conduc- 
tivity through k — vlC/3, being v the Fermi velocity of 
the electrons and I their mean free patbJ^. We plot this 
quantity in the right upper panel of Fig. [3) This phys- 
ical quantity is almost insensitive to the opening of the 
energy gap and the different plots, corresponding to dif- 
ferent values of At almost coincide within the scale of 
the figure. From the lower panel of Fig. [3] we see that 
while for high temperatures there is a linear relation be- 
tween K and C , this is not the case at lower temperatures 
where Andreev type processes are relevant. 

As stressed before, our calculation is not restricted to 



small temperature gradients. We show in Fig. [4] a plot 
of the heat current for several values of the anisotropy 
parameter At as a function of the temperature of the 
left chain while the temperature of the right chain T is 
set fixed to zero. The figure clearly shows the suppression 
of the current as a consequence of the "Andreev reflection 
"phenomena mentioned before. In fact for T < Al the 
current is exponentially small, while it grows for higher 
temperatures. 

Finally, in Fig. \5\ we illustrate the behavior of the 
heat current when finite different magnetic fields are ap- 
plied at the two side chains. The effect of applying mag- 
netic fields at both sides of the junction leads to a in- 
teresting effect which me name "thermal diode effect". 
As discussed in the previous section, a finite magnetic 
field originates a shift in the arguments of the functions 
(w), which leads to asymmetries in the transmission 
functions T n (cu) and T a (cu). For At = 0, only the normal 
transmission function and the functions Y^ 9 (lo) are non- 
vanishing. Furthermore, these functions are identical and 
gapless for a — L,R. Therefore, the heat flow is per- 
fectly antisymmetrical (the sign of the current is reversed 
preserving the absolute value) under the simultaneous 
change Tl <-> Tr and Bl — * Br. Instead, for a finite 
Al, the situation changes. A gap opens for the excita- 
tions of the left chain and the functions Y V L V (u>) = van- 
ish for \u\ < Al, while the functions T R 3 (u) = Y R (—u) 
remain finite. The consequence is an asymmetry in the 
behavior of the transmission functions under the change 
Bl — ► Br. The result is an effect of thermal rectifica- 
tion. That is, the magnitude of the current J® when 
(Tl = T,ial = jtx) and (Tr = T 1 ,\ir — //) is different 
to J Q when (T L = T',fi L = and (Tr = T,fj, R = n), 
which means that the device is more likely to conduct 
heat when the temperature difference is applied in one 
direction than in the other. We display the phenomena 
for two different values of A t . We show the current when 
Oi = 0.3, Tl = T,fj, R = 0,T R = 0) with dots and the 
current when (/j,r = 0.3, Tr =T,/xl = 0, Tl = 0) with a 
full line. When the value of At = 0.2 both currents are 
rather large and similar but when Al = 0.75 the currents 
are smaller and clearly different. 

IV. SUMMARY AND DISCUSSION 

We have presented a theoretical framework to study 
heat transport in one-dimensional spin heterostructures. 

In the present work we have focussed on a simple sys- 
tem composed of a junction between an anisotropic (XY) 
and an isotropic (XX) chain under the effect of an inho- 
mogeneous magnetic field along the z direction. Using 
the Jordan- Wigner transformation to map the problem 
into a fermionic system and using the non equilibrium 
Keldysh-Schwinger formalism we have obtained exact ex- 
pression for the heat current in terms of Green's functions 
of the " disconnected" spin chain components. The result- 
ing expressions can be evaluated numerically in a simple 
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some of these issues in future work. 
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APPENDIX A: GREEN'S FUNCTIONS FOR AN 
OPEN CHAIN WITH p-WAVE 
SUPERCONDUCTIVITY 



FIG. 5: (Color online) Thermal Diode effect. We show the 
heat current when (fiL = 0.3, Tl = T, fj,R = 0, Tr = 0) (dot- 
ted line) and when (pt R = 0.3, Tr = T, fi L = 0,T L = 0) for 
two values of the pairing parameter Al = 0.2, 0.75 



way. In the limits A L < T < J and T < A L < J 
explicit analytic expressions can also be given. We have 
studied the heat transport as a function of the different 
parameters of the model and we have shown that when 
different magnetic fields are applied at the end chains, a 
rectifying effect in the heat current occurs. This effect 
might be of interest for applications. Its origin can be 
traced back to the appearance of paring terms induced 
by the anisotropy parameter, which are in turn responsi- 
ble for an Andreev reflection type mechanism. 

In this work we have analyzed a simple model. How- 
ever this methodology can be straightforwardly extended 
to more complex structures with many junctions and dis- 
order. Our treatment relies on the Jordan- Wigner trans- 
formation which maps the original spin Hamiltonians into 
fermionic ones. In the case we have considered, the lat- 
ter are bilinear. In more generic models, although we 
expect the rectifying effect still to be present, the tech- 
nical analysis could be more complicated. For instance, 
in the isotropic Heisenberg model, which in addition to 
the exchange interaction along x and y directions, con- 
tains an additional exchange term along the z direction, 
the Jordan- Wigner transformation translates such a term 
into a many-body fermionic interaction, which does not 
enable a straightforward analytical solution of the prob- 
lem, as in the case we considered here. Nevertheless, the 
Green's function formalism offers a framework for the 
construction of systematic approximations to treat those 
terms. Numerical methods could also be useful to deal 
with models containing many-body terms ^ As in elec- 
tronic systems, many-body terms are expected to intro- 
duce further inelastic scattering processes, which could 
add further ingredients in addition to the transport mech- 
anisms we have discussed here. We hope to report on 



In this appendix we show a derivation of the Green's 
functions G L11 (u), G l11 (uj), F£ M (w) and F l11 {uj), 
entering Eq. (|27[). which correspond to the end of a 
half infinite chain with p-wave superconductivity in the 
BCS approximation. Making the superconducting pa- 
rameter Al — 0, the first two Green's functions give the 
corresponding result for the normal chain g® R l and 

~g\i i(u>) respectively. 

The Green's functions of the open chain can be solved 
considering a ring of N sites, periodic except for the fact 
that the energy at one site (which we label as site 0) 
is increased by an energy A, and then taking the limit 
N, A — ► +oo . The Hamiltonian is 

N-l 

h = £{M///*+i + //+i/*] 

1=0 

+A[//// +1 +/ m / i ]} 

-V E + A /o/o- (Al) 

1=0 

We have solved the problem using two different meth- 
ods: i) solving the equations of motion in Fourier space, 
and ii) solving a Dyson's equation that relates the above 
Green's functions to those of the periodic chain (A = 0) 
which can be obtained easily using Bloch theorem. Both 
results of course coincide, but the latter method involves 
a simpler algebra. We define a matrix 

- / G? n (uj) F? , , (w) \ 
G=( J' 1 ' 1 ) J' 1 ' 1 ) , (A2) 

with the Green's functions for A ^ 0, and a correspond- 
ing matrix g for A = 0. These matrices are equivalent 
to the ones obtained by using Nambu's representation 
for the Hamiltonian and the Green's functions^ From 
the equations of motion of these Green's functions, one 
obtains 
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G = g + gAG, 

where A, is proportional to A. 

Solving Eq. ((A3)) for G and taking the limit A 
the following expressions result 



h (uj) 
h 2 (oj) 



h (uj) h* (-uj) 
ho(-u>) h (uj) 



(A3) 



(A4) 
(A5) 



For N — > +00, the sums can be replaced by integrals. 
Decomposing the integrands into a sum of simple frac- 
tions with denominators linear in cos k and numerators 
independent of fc, the integrals can be evaluated analyt- 
ically using^i 



n Jo 



dk 



cosk + b V6 2 - 1 ' 



(A8) 



where the sign of the root is determined by the sign of 
the imaginary part of the second member. 
Defining 



The h- functions entering the second members of Eqs. 
(|A4|) and (|A5jl are Green's functions of the periodic chain 
and can be calculated easily in Fourier space. The result 
is 



h Q (uj) 
hi(u>) 
h 2 (oj) 



-y- 

k 

1 E 

k 



uJ + e k 



(u + Cjfe) cosfc 



2A sin 2 k 



N 



(A7) 



where in the second members to includes an infinites- 
imally small imaginary part, e k = 2wcosk — fj,, and 
A fc = 2Asinfc. 



u z A _ a . X9 
uj = — , A = — , (U = — , d=l — A , 

r = [(£ 2 - A 2 )d+(A/~i) 2 j 1/2 /rf- 
61 = jl/d — r, b 2 = fl/d+ r, 
di = [4wrd] 1 , 

the result takes the form 



(A9) 



h (uj) = di[{r - & - fth 2 /b)I{bi) 

+ (r + u + fiA 2 /b)I(b 2 )}, 
hi(u) = -d x {2r- [(r-w-/iA 2 /6)6i/(6i) 

+ (r + ^ + /iA 2 /6)6 2 /(6 2 )]}, 
fta(w) = diA{2r+[/(6 1 )]- 1 -[/(6 2 )]- 1 }. (A10) 
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